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A three-dimensional layered Ising-Antiferromagnet with a ferromagnetic intra-layer 
coupling to 2 neighbors, zj > 0, and an antiferromagnetic interlayer coupling to 
z' neighbors, z' J' < 0, is investigated by Monte Carlo simulations on a hexag- 
onal lattice. The physical nature of the anomalous temperature bahavior of the 
sublattice magnetizations, which is found for certain values of r = zJ/z'J' and 
2' in magnetic fields is explained in terms of successive phase transitions. They 
take place on the ferromagnetic 2-dimensional spin-down sublattice at T ss T)?'*, 
smeared by a finite stabilizing molecular field, and on both antiferromagnetically 
coupled sublattices at T^'' > T^'*. 



1 Introduction 

Anisotropic three-dimensional antiferromagnets (AFs) are still an interesting 
subject of theoretical investigations. There exist many investigations e.g. on 
their magnetic pha|aa_diagram within the framework of Ising or anisotropic 
Heisenberg models.ErlHl It is found, that layered Ising AFs with two compet- 
ing interaction parameters (Fig. |l|) may exhibit two rather different phase 
diagrams (Fig. g) depending only on two parameters, the ratio r ~ zJ/z'J' 
and z'Jj where J and J' are the coupling constants of the intra-layer and of 
the inter-layer exchange, respectively; z and z' are the coordination numbers 
of the couplings. The Ising Hamiltonian is of the form: 

n=^j J2 ^^^1 -■^' Yl ^^^^ -hJ2s,, (1) 

where H is the applied magnetic field acting on all spins Si, with Si = ±1. The 
two kinds of phase diagrams are continuously transformed into one another 
by changing the crucial parameters r and z!. For large values of |r|, i^ 
Irl > 0.6J3 and small values for z', i.e. z' < 10,a the FeCl2-like phase diagramtj 
is found (Fig. (b)), whereas for low values of \r\ and large values for z' 
the case (b) in Fig. appears. This second case is characterized by three 
interesting features. On the one hand, a possible decoupling of the tricritical 
point (TCP) into a critical endpoint (CEP) and a bicritical endpoint (BCE) 
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Figure 1. Schematic drawing of a layered 3-dimensionaI antiferromagnet with a ferromag- 
netic coupling in the layers, Ja—A = Jb—b = zj , and an antiferromagnetic coupling 
between adjacent layers, Ja—b = Jb—A = z'J'. The magnetic field, H, is applied along 
the antiferromagnetic stacking direction. 



is encountered. From previous investigatioHS it follows, that this decoupling is 
only observed in mean field calculati(ms,Q"Ll while in Monte Carlo simulations 
only one multicritical point is foundJlj On the other hand, the second-order 
phase line has a balloon-like shape and extends even above the limiting field 
value Hco (spin-flip field at T = 0). This means, that for some fixed field 
values H > Hco it is possible to cross this phase line twice with increasing 
temperature T. Furthermore, above and below the critical line anomaly lines 
are found, where the magnetization exhibits an additionaLinflection point and 
the specific heat shows an additional broad maximum. BEj 
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Figure 2. Schematic magnetic phase diagrams of 3-dimensional layered Ising AFs (see Fig- 
ure hi) choosing (a) \r\ large and z' small (e.g. \r\ = 1.0 and z' = 4), and (b) vice versa (e.g. 
\r\ = 0.4 and z' = 20). Bold, thin and broken lines refer to first- and second-order phase 
transitions and anomalies, respectively. Transitions occur between the antiferromagnetic 
(AF) and paramagnetic (PM) phase. TCP, CEP, BCE denote tricritical, critical end- and 
bicritical endpoints, respectively. 
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Figure 3. (a) H-T phase diagram of a 3-dimensionaI Ising antiferromagnet on a hexagonal 
lattice (24x24x18) with coupling constants zJ = 6 ■ 0.7 and z' J' = -20 • 0.5 (fcs = 1). Pe- 
riodic boundary conditions for all axes were used. For T < Tmcp ~ 3 the phase transition 
is discontinuous, while for T > Tmcp ^ second-order phase transition occurs at the line 
Hc{T). The broken lines denoted as H— and H^ are the anomaly lines, (b) Calculated 
magnetization M (open circles) and its derivative dM/dT (full circles) vs temperature T in 
a field H = 9.5. 



2 Comparison between theory and experiment 

Fig. II (a) shows the H-T phase diagram of a 3-dimensional layered Ising AF 
obtained from Monte-Carlo simulations on a hexagonal lattice with periodic 
boundary conditions. All quantities like T, H, J and J' are considered to 
be dimensionless. For the simulation the Metropolis algorithm was used with 
/cb = 1. The parameters z, z', J, J' are chosen such as to reproduce case (b) of 
Fig. 0. One observes only one multicritical point (MCP), where a first-order 
phase line (T < Tmcp) and a second order phase line, Hc{T) meet, and where 
the two anomaly lines, iJ_ and iJ+, originate. 

Fig. 13 (b) shows one magnetization curve, referring to the phase diagram 
of Fig. H^(a), as a function of temperature for H = 9.5 and its derivative, 
dM/dT. The magnetization is defined to be M ~ ^i^i/N, where N is 
the number of lattice sites. M vs T clearly shows an anomalous curvature 
for T < Tc, which manifests itself in the derivative as an additional broad 
maximum. 

The phase diagram shown in Fig. y (a) resembles^JJiat 
sulating uniaxial antiferromagnet with Tn = 14.2 K Ej E^ 



of FeBr2, an in- 
(Fig. I). Both 
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Figure 4. (a) H-T phase diagram of FeBr2, where H is the apphed axial magnetic field. 
The lines Ha,Hci,Hi^2,H— and i?+ were obtained by magnetization or susceptibility mea- 
surements. H\ and H'^ denote the positions of the spikes found by specific heat (open 
diamonds) (Ref. hj) and off-axis magnetization measurements vs temperature (solid cir- 
cles), respectively, (b) Magnetization M vs temperature T and its derivative for a field 
H = 2.07 MA/m (field parallel to the c-axis). The phase lines H\ or H'^ are not seen in this 
configuration (see Ref. hq). The broad peak at T_ = T{H-^) is due to strong non-critical 
fluctuations (Ref. hsl). 



in Fig. 1^ and Fig. |4| lines of non-critical fluctuations (or anomaly lines) do 
appear. In this article the attention is focused on these non-critical fluctua- 
tions at T_ p.ll(|££_). Although this phenomeapn is well investigated both in 
experimental3't3ll3 on FeBr2 and in theoryEillJ no clear explanation for the 
occurrence of these fluctuations yet exists. Especially one wonders, how it is 
possible that two completely different types of phase diagrams are found by 
varying only two parameters, r and z' . 

3 Non-critical fluctuations 

In order to gain insight into the origin of the non-critical fluctuations at 
H-{T), we performed systematic Monte Carlo simulations of the sublat- 
tice magnetizations with the same parameters as in Fig. H, zJ = 4.2 and 
z'J' — —10.0. The exchange constants and especially the number of cou- 
pled neighbours are comparable to those in FeBr2.EZI However, here we used 
only two exchange couplings. The different intra-planar exchange parameters 
found in FeBr2 were absorbed in one ferromagnetic effective intra-planar cou- 
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Figure 5. Sublattice magnetization of an Ising AF with an hexagonal lattice (24x24x18) and 
zj = 6 0.7 and z' J' = —20 ■ 0.5 (kg = 1). For low values of H the sublattice magnetization 
Ma and Mb are as expected for an uniaxial antiferromagnet. By increasing the field value 
the sublattice magnetizations become strongly asymmetric. 



pling constant. Fig. H shows the temperature dependences of the sublattice 
magnetizations, Ma and AIb (see Fig. ^, for H — 0, A, 8, 9.5 and 9.95, 
respectively. While Ma (parallel) and Mb (antiparallel to H) are symmetric 
for zero field, Ma{T) — —Mb(T), they become more and more inequivalent 
with increasing field. For fields coming close to the spin-flip one, Hcq — fO, 
anomalous bumps appear at T « T- in Mb{T), whereas Ma{T) is virtually 
constant up to that temperature. Obviously all of the non-critical fluctuations 
observed at T « T- happen to occur merely on the B-sublattice. In other 
words, these fluctuations are essentially constrained to 2-dimensional (2d) lay- 
ers separated by magnetically saturated up-spin layers of the A-sublattice. It 
is, hence, tempting to compare the B-sublattice with an ensemble of 2d ferro- 
magnets (FMs) with the same intra-layer parameters, zJ > 0, but subjected 
to a field Heff = H — H^o, where H is the field applied to the corresponding 
3d AF. 

Fig. @ shows the sublattice magnetization curves for H — 9.5 and 9.95 
as before and the magnetization curves of a 2d FM with zJ — 4.2 in fields 
Heff — —0.5 and —0.05, respectively. One observes, that the magnetization 
curve of the 2d FM fits well with the spin-down sublattice magnetization Mb 
up to the inflection point at T_ (arrows) . It seems, that the magnetization 
of the spin-down sublattice behaves like a 2d FM in an effective mean field. 
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Figure 6. Two examples of sublattice magnetizations (circles) of the model as in Fig 
{HcO = —z' J' = 10.0) and in addition magnetization curves of a 2d FM (solid lines) wit 



zj = 6- 0.7 in fields H, 



eff 



H — Hrj 



0, —0.05, and —0.5, respectively. Arrows denote 



the inflection point, which is associated with T_ 



which is the sum of the external applied field H and the field, produced by 
the fully magnetized spin- up sublattice (A-sublattice), H^ff = H + Ha- Since 
Ma is completely magnetized, Ha = z'J' = — i?cO, one has Heff = H — Hcq- 
If the effective field H^ff becomes zero {Heff = '^ H — Hcq) we have 
the case of a 2d FM in zero field, which undergoes a phase transition at 
T^'^ and becomes paramagnetic for T > T^"^ (see Fig. ||). Although this 
case cannot be obtained here, because of the spin-flip occurring in the range 
Hmcp < H < HcO, the anomaly temperatures T_ can be associated with the 
points of inflection of the AI{Heff) vs T curves of the 2d Ising FM (arrows 
in Fig. g). Therefore one can conclude, that the anomaly of Mb singnifies 
the thermal destruction of 2d ferromagnetic order on the quasi-decoupled 
B-sublattice layers, which precedes the global 3d phase transition of both 
sublattices. In other words, the anomaly we find in the magnetization is due 
to the finite temperature range lying between T^"^ and T^'^. If this splitting 
is reduced by increasing the intra-planar ferromagnetic interaction and thus 
increasing the 2d transition temperature, the anomaly decreases and vice 
versa. 
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4 Conclusion 

The occurrence of anomalies, which are observed in certain 3d Ising AFs is 
due to the separation of the smeared 2d phase transition on one sublattice 
from the 3d global phase transition of both sublattices. This is typical of an- 
tiferromagnets, whose inter-planar antiferromagnetic coupling, is strong com- 
pared with the intra-planar exchange {\r\ — \zJ/z'J'\ < 0.6)13. A high value 
of the inter-planar exchange gives rise to a high value of the spin-flip field 
HcQ = —z'J'. This results in a strong stabilization of the spin-up sublattice, 
which acts only as a mean-field Hmf ~ z' J' on the spin-down sublattice. 
Hence, the spin-down sublattice behaves like a 2d FM in an effective field 
Heff — H + Hmf ■ The quality of the mean- field approximation crucially de- 
pends on a high coordination number z' , which is needed for integrating out 
thermal spin-fluctuations. 
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